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[i]  During  spring  and  summer  time,  coastal  upwelling  influences  circulation  and  ecosystem 
dynamics  of  the  Monterey  Bay,  California,  which  is  recognized  as  a  National  Marine 
Sanctuary.  Observations  of  physical,  bio-optical  properties  (including  bioluminescence) 
together  with  results  from  dynamical  biochemical  and  bioluminescence  models  are  used 
to  interpret  the  development  of  the  upwelling  event  during  August  2003  in  Monterey 
Bay,  California.  Observations  and  the  biochemical  model  show  the  development  of  a 
phytoplankton  bloom  in  the  southern  portion  of  Monterey  Bay.  Model  results  show  an 
increase  of  nutrients  in  the  southern  portion  of  the  bay,  where  nutrient-rich  water  masses  are 
brought  in  by  the  southward  flow  and  cyclonic  circulation  inside  the  bay.  This  increase 
in  nutrients  together  with  the  sluggish  circulation  in  the  southern  portion  of  the  bay  provides 
favorable  conditions  for  phytoplankton  growth.  Our  observations  and  models  suggest 
that  with  the  development  of  upwelling  the  offshore  water  masses  with  the  subsurface 
layer  of  bioluminescent  zooplankton  were  replaced  by  water  masses  advected  from 
the  northern  coast  of  the  bay  with  a  relatively  high  presence  of  mostly  nonbioluminescent 
phytoplankton.  Inshore  observations  from  autonomous  underwater  vehicles  (AUVs)  show 
consistent  coincidence  of  chlorophyll,  backscatter,  and  bioluminescence  maxima  during 
upwelling  development.  Offshore  AUV  observations  (taken  at  the  entrance  to  the  bay) 
show  a  deeper  bioluminescence  maximum  below  the  surface  layers  of  high  chlorophyll 
and  backscatter  values  during  the  earlier  stages  of  upwelling  development.  Later,  the 
observed  deep  offshore  bioluminescence  maximum  disappeared  and  became  a  shallower 
and  much  weaker  signal,  coinciding  with  high  chlorophyll  and  backscatter  values  offshore. 

Based  on  the  biochemical  and  bioluminescence  models,  a  methodology  for  estimating 
the  nighttime  water-leaving  radiance  due  to  stimulated  bioluminescence  is  demonstrated 
and  evaluated. 

Citation:  Shulman,  1.,  M.  A.  Moline,  B.  Penta,  S.  Anderson,  M.  Oliver,  and  S.  H.  D.  Haddock  (201 1),  Observed  and  modeled 
bio-optical,  bioluminescent,  and  physical  properties  during  a  coastal  upwelling  event  in  Monterey  Bay,  California,  J.  Geophys. 
Res .,  116 ,  C01018,  doi:  10.1 029/20 10JC006525. 


1.  Introduction 

[2]  A  multi-institution,  multidisciplinary  Autonomous  Ocean 
Sampling  Network  field  experiment  (called  AOSN  II)  was 
conducted  in  the  Monterey  Bay,  California,  during  August 
and  September  2003  [Ramp  et  al. ,  2009].  The  Monterey  Bay 
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represented  an  attractive  site  for  the  experiment  due  to  the 
availability  of  an  already  existing  observing  system  including 
the  established  network  of  HF  radars,  semipermanent  moored 
surface  buoys  and  communications  infrastructure.  During 
northwesterly,  upwelling  favorable  winds,  the  mesoscale 
and  large-scale  features  in  and  around  the  Monterey  Bay 
are  mostly  determined  by  the  interaction  between  upwelling 
filaments  formed  at  headlands  to  the  north  (Pt.  Ano  Nuevo, 
Figure  1 )  and  to  the  south  of  the  Bay  (Pt.  Sur)  and  by  the 
California  Current  system  offshore  of  the  Bay  [Rosenfeld 
et  al ,  1994],  Wind-driven  upwelling  of  nutrient-rich  water 
supports  high  levels  of  phytoplankton  productivity  in  Mon¬ 
terey  Bay  [Pennington  and  Chavez ,  2000;  Ryan  et  al,  2005]. 
Late  stage  phytoplankton  blooms  yield  an  accumulation 
of  autotrophic  dinoflagellates  that  include  bioluminescent 
species  [Haddock  et  al,  2010;  Moline  et  al,  2009]. 
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Figure  1.  (a)  The  NCOM  ICON  model  domain  (bounded  by  solid  black  line)  and  Ml  mooring  location, 
(b)  Locations  of  HF  radar  sites;  V-shaped  transect  of  CalPoly  AUV  REMUS  and  sections  sampled  by 
AUV  DORADO,  (c)  Observed  wind  velocities  at  Ml  during  August  2003. 


[3]  One  objective  of  the  2003  experiment  was  to  learn  how 
to  apply  new  tools,  technologies,  and  analysis  techniques 
to  adaptively  sample  the  coastal  ocean.  Another  objective 
was  to  use  the  gathered  information  to  develop  accurate 
forecasts  of  the  bay-scale  patterns  of  physical  and  biological 
fields,  including  bioluminescence  [Ramp  et  al. ,  2009; 
Haddock  et  al. ,  2010;  Moline  et  al. ,  2005,  2009;  Shulman 
et  al .,  2005].  For  this  reason,  the  field  program  included  an 
extensive  sampling  of  the  bay  and  surrounding  areas  with  a 
fleet  of  underwater  gliders,  propeller-driven  autonomous 
underwater  vehicles  (AUVs),  a  low-flying  aircraft  and  HF 
radars  in  addition  to  moorings,  ships,  and  other  more  tradi¬ 
tional  observational  techniques  (see  also  special  issue  of 
Deep  Sea  Research  Part  //,  56 ,  pp.  61-260,  2009). 

[4]  Modeling  of  physical  conditions  during  the  experiment 
was  presented  by  Shulman  et  al.  [2009,  2010].  The  model 
(based  on  the  Navy  Coastal  Ocean  Model  (NCOM))  was  able 
to  reproduce  many  observed  circulation  features  of  upwell- 
ing  and  relaxation  events  of  August  2003.  The  physical 


model  is  coupled  to  the  biochemical  submodel  of  Chat 
et  al.  [2002].  Predictions  from  the  physical  model  are  also 
used  to  simulate  changes  in  the  bioluminescence  intensity  in 
accordance  with  methodology  outlined  by  Shulman  et  al. 
[2003,  2005]. 

[5]  The  primary  objective  of  the  present  study  is  to 
describe  and  provide  the  interpretation  of  the  upwelling 
development  during  the  experiment.  For  this  reason,  obser¬ 
vations  of  physical,  bio-optical  properties  (including  biolu¬ 
minescence)  together  with  results  from  physical,  biochemical 
and  bioluminescence  models  are  used  to  achieve  this  goal. 
Also,  biochemical  and  bioluminescence  dynamical  models 
are  combined  together  to  demonstrate  predictions  of  the 
nighttime  water-leaving  radiance  due  to  stimulated  subsur¬ 
face  bioluminescence. 

[6]  The  structure  of  the  paper  is  as  follows:  Section  2 
describes  observations  and  models.  Analysis  of  observed 
and  modeled  bio-optical,  physical  properties  during  the 
upwelling  event  is  presented  in  section  3.  Section  4  is  devoted 
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to  the  bioluminescence  and  the  nighttime  water-leaving  radi¬ 
ance  modeling.  Conclusions  are  presented  in  section  5. 

2.  Methods 
2.1.  Observations 

[7]  Observations  of  winds  and  water  velocity  from  the 
Monterey  Bay  Aquarium  Research  Institute  (MBARI)  sur¬ 
face  mooring  Ml  (36.74°N,  122.02°W)  are  used  in  this 
study  (Figure  1). 

[s]  During  the  2003  field  program,  the  mooring  had  a 
downward  looking  RD  Instruments  Inc.  75  KHz  Acoustic 
Doppler  Current  Profiler  (ADCP)  which  was  set  up  to 
sample  currents  every  15  min.  In  60  8  m  bins  up  to  500  m 
depth  (the  first  bin  at  16  m  depth).  Surface  wind  speed  and 
direction  were  measured  by  a  RM  Young  model  05103 
wind  monitor. 

[9]  Surface  current  observations  used  in  this  study  were 
derived  from  a  network  of  SeaSonde-type  HF  radar  instru¬ 
ments  deployed  in  the  Monterey  Bay  region.  These  instru¬ 
ments  exploit  information  in  the  radio  wave  backscatter 
from  the  ocean  surface  to  infer  movement  of  the  near  surface 
water.  Each  individual  SeaSonde  instrument  provides  a 
distribution  of  “radial”  velocity  observations  each  hour  on  a 
polar  coordinate  grid  centered  on  the  radar  site.  Vector 
currents  were  estimated  on  a  Cartesian  grid  with  a  horizontal 
resolution  of  3  km  by  computing  the  best  fit  vector  velocity 
components  using  all  radial  velocity  observations  within  a 
radius  of  3  km  for  each  grid  point  each  hour  [Paduan  and 
Shulman ,  2004].  During  the  AOSN  II  experiment,  surface 
currents  were  estimated  based  on  input  from  four  HF  radar 
sites  (Figure  1):  Santa  Cruz  (SCRZ),  Moss  Landing 
(MLNG),  The  Naval  Postgraduate  School  (NPGS)  and  the 
Point  Pinos  (PPIN). 

[10]  The  goal  of  AUVs  surveys  was  sampling  of  biolu- 
minesccnce  (BL).  BL  was  measured  using  two  MBBP 
series  custom  bathyphotometers  developed  at  University  of 
California  Santa  Barbara.  Bathyphotometers  were  mounted 
on  a  CalPoly  AUV  (Remote  Environmental  Monitoring 
Units  (REMUS))  [Moline  et  al. ,  2005]  and  a  larger  MBARI 
DORADO  AUV  [Haddock  et  al .,  2010].  The  biolumines- 
ccnce  bath ypho tom ctcr  [Herren  et  al,  2005]  pumps  water 
into  a  0.5  L  sample  chamber  at  a  rate  of  0.35-0.4  L/s.  Flow 
rate,  temperature,  and  light  levels  (photon  flux,  assuming 
isotropic  emission)  are  measured  in  the  sample  chamber.  The 
instrument  is  calibrated  both  radiometrically  by  a  known 
light  source,  and  biologically  by  insertion  of  a  known  con¬ 
centration  of  dino flagellates.  With  the  bathyphotometer,  we 
capture  a  repeatable  and  representative  fraction  of  the  bio¬ 
luminescence  present  in  the  environment.  The  CalPoly  AUV 
conducted  surveys  along  a  V-shaped  transect  in  the  northern 
half  of  the  Monterey  Bay  (Figure  lb).  These  operations 
began  near  Santa  Cruz,  ran  out  to  MBARI  surface  buoy  Ml , 
then  returned  back  to  shore.  These  REMUS  runs  collected 
data  between  2 1 00  and  0400  LT  each  night  on  10-18  August 
2003.  The  vehicle  followed  a  sawtooth  pattern  to  40  m  depth 
sampling  with  a  CTD,  transmisso meter,  fluorometer,  and 
ADCP  in  addition  to  a  bio  luminescence.  The  DORADO  AUV 
sections  are  shown  on  Figure  I .  Instruments  on  board  included 
a  CTD,  fluorometer,  oxygen  and  nitrate  sensors,  biolumines¬ 
cence,  and  ADCP. 


2.2.  Models 

[11]  The  Monterey  Bay  model  (called  the  NCOM  ICON) 
consists  of  the  physical  model  [Shulman  et  al .,  2007],  which 
is  coupled  to  the  biochemical  model  [Chai  et  al .,  2002]. 
The  physical  model  of  the  Monterey  Bay  is  based  on  the 
NCOM  model,  which  is  a  primitive  equation,  3-D,  hydro¬ 
static  model.  It  uses  the  Mellor-Yamada  level  2.5  turbulence 
closure  scheme,  and  the  Smagorinsky  formulation  for  hori¬ 
zontal  mixing  [ Martin ,  2000].  The  biochemical  model  of  the 
NCOM  ICON  simulates  dynamics  of  two  sizes  of  phyto¬ 
plankton,  small  phytoplankton  cells  (<5  //m  in  diameter)  and 
diatoms,  two  zooplankton  grazers,  nitrate,  silicate,  ammo¬ 
nium,  and  two  detritus  pools  [Chai  et  a l .,  2002]. 

[12]  The  NCOM  ICON  model  is  set  up  on  a  curvilinear 
orthogonal  grid  with  resolution  ranging  from  I  to  4  km.  The 
model  domain  is  shown  on  Figure  1 .  The  model  is  forced  with 
surface  fluxes  from  the  Coupled  Ocean  and  Atmospheric 
Mesoscale  Prediction  System  (COAMPS)  [Doyle  et  al .,  2009] 
at  3  km  horizontal  resolution.  The  3  km  resolution  CO  AM  PS 
grid  mesh  is  centered  over  Central  California  and  the  Mon¬ 
terey  Bay.  Phytoplankton  photosynthesis  in  the  biochemical 
model  is  driven  by  Photosynthetically  Active  Radiation 
(PAR),  which  is  estimated  based  on  the  shortwave  radiation 
flux  from  the  COAMPS  model.  The  Penta  et  al.  [2008] 
scheme  is  used  for  PAR  attenuation  with  depth. 

[13]  The  NCOM  ICON  model  uses  the  Navy  Coupled 
Ocean  Data  Assimilation  (NCODA)  system  [Cummings, 
2005]  for  the  assimilation  of  the  temperature  and  salinity 
data  from  different  observational  platforms.  The  NCODA  is  a 
fully  3-D  multivariate  optimum  interpolation  system.  Assi¬ 
milation  of  temperature  and  salinity  data  is  performed  every 
12  h  (assimilation  cycle).  Differences  between  the  NCODA 
analysis  and  the  model  forecast  are  uniformly  added  to  the 
model  temperature  and  salinity  fields  over  the  assimilation 
cycle.  The  evaluation  of  the  NCOM  ICON  model  predictions 
are  presented  by  Shulman  et  al.  [2007,  2009,  2010]. 

[14]  Open  boundary  conditions  for  the  NCOM  ICON  are 
derived  from  the  regional  model  of  the  California  Current 
(NCOM  CCS)  [Shulman  et  al.,  2007].  The  NCOM  CCS  has  a 
horizontal  resolution  of  about  9  km  and,  the  model  is  forced 
with  atmospheric  products  derived  from  the  CO  AMPS 
[Doyle  et  al.,  2009]. 

[15]  Open  boundary  conditions  for  the  regional  NCOM 
CCS  model  are  derived  from  the  NCOM  global  model 
[Rhodes  et  al,  2002;  Barron  et  al,  2004],  which  has  1/8° 
horizontal  resolution.  The  model  assimilates  satellite-derived 
sea  surface  height  (SSH)  and  sea  surface  temperature  (SST) 
data  via  synthetic  temperature  and  salinity  profiles  derived 
from  the  Modular  Ocean  Data  Assimilation  System  (MODAS) 
[Foxetal,  2002],  and  uses  atmospheric  forcing  from  the  Navy 
Global  Atmospheric  Prediction  System  (NOGAPS)  [Rosmond 
et  al. ,  2002]. 

[16]  The  bioluminescence  model  (BL  model)  is  based  on 
BL  predictions  with  an  advection-diffusion-reaction  model 
(ADR),  with  velocities  and  diffusivities  taken  from  the 
NCOM  ICON  model  [Shulman  et  al,  2005,  2003].  The 
BL  model  consists  of  dynamical  initialization  of  the  ADR 
model  by  assimilating  available  BL  observations,  and  fore¬ 
casting  BL  intensity  with  the  ADR  model.  There  is  no 
modeling  of  sources  and  sinks  of  BL  intensity  due  to  bio¬ 
logical  interactions,  nor  modeling  of  behavioral  dynamics 


3  of  13 


CO 10 18 


SHULMAN  ET  AL.:  BIO-OPTICS,  PHYSICS  DURING  UPWELLING 


C01018 


of  bioluminescent  organisms  (for  example  vertical  migration 
of  dinoflagellates). 

3.  Bio-Optical  and  Physical  Properties  During 
the  Upwelling  Event 

[17]  In  accordance  with  the  observed  wind  velocity  at 
the  mooring  Ml,  AUVs  REMUS  and  DORADO  surveys 
were  conducted  during  the  extended  upwelling  event  of 
7-19  August  (Figure  1).  The  spatial  distribution  of  the  HF 
radar  derived  surface  currents  and  the  subsurface  profiles  of 
northward  and  eastward  velocity  components  at  mooring  Ml 
are  shown  on  Figure  2.  Both  surface  and  subsurface  currents 
are  averaged  over  3  days  of  upwelling  (15-17  August). 
Figure  2  indicates  that  during  the  upwelling  event,  there  was 
a  development  of  strong,  wide  southward  flow  at  the  surface, 
which  extends  up  to  1 50m  depth  (Figure  2,  negative  values  of 
V  component  of  the  mooring  M 1  currents  indicate  southwfard 
flow).  This  southward  flow  separates  a  pair  of  cyclonic 
(inside  the  bay)  and  anticyclonic  (outside  the  bay)  circula¬ 
tions.  The  satellite-derived  SST  and  chlorophyll  images 
(from  MODIS-AQUA  satellite)  indicate  that  southward  flow 
along  the  entrance  to  the  bay  has  colder  surface  temperatures 
than  water  masses  offshore  and  inside  the  bay  (Figure  2). 
Also,  there  is  a  lower  surface  chlorophyll  concentration  in  the 
southward  flow  than  that  observed  inside  the  bay.  The  off¬ 
shore  anticyclonic  circulation  has  warmer  and  less  productive 
surface  water  masses  in  the  center  of  the  circulation  and 
relatively  more  productive  water  masses  at  the  eastern  flank 
of  the  circulation  (western  side  of  the  southward  jet).  In  the 
bay,  there  are  surface  water  masses  with  warmer  and  higher 
chlorophyll  values  along  the  coast.  Satellite  images  indicate 
the  development  of  a  phytoplankton  bloom  in  the  southern 
part  of  the  bay  which  extends  from  Moss  Landing,  California 
toward  the  south. 

[is]  REMUS  surveys  on  11  and  12  August  (Figure  3) 
show  coincidence  of  inshore  maxima  of  chlorophyll,  BL  and 
backscatter  distributions.  This  suggests  that  the  inshore  BL 
maximum  is  associated  with  the  planktonic  community 
(dinoflagellates).  This  conclusion  was  supported  by  using 
general  differences  in  flash  kinetics  between  planktonic 
dinoflagellates  and  zooplankton  in  Moline  et  al.  [2009].  The 
REMUS  observations  also  show  a  strong  BL  signal  offshore 
in  the  deeper  water  around  the  Ml  mooring  location  (center 
of  the  entrance  to  the  bay)  (Figure  3).  This  deeper  offshore 
BL  maximum  is  below'  the  observed  chlorophyll  layer. 
Low  values  in  optical  backscatter  in  the  area  of  the  BL 
maximum  (Figure  3)  suggest  that  this  observed  deeper  BL 
signal  is  due  to  larger  zooplankton  and  probably  not  due  to 
heterotrophic  dinoflagellates. 

[19]  Strong  correlations  between  inshore  maxima  in  chlo¬ 
rophyll,  BL  and  backscatter  persists  over  the  next  days  of 
REMUS  sampling  (13-15  August)  (Figure  3)  and  in  the 
DORADO  survey  on  1 3  August  (going  from  the  M 1  mooring 
to  the  south  of  the  bay,  Figure  4).  At  the  same  time,  the  off¬ 
shore  deep  BL  signal  almost  disappeared  on  1 3  August  in  area 
around  mooring  Ml  in  the  REMUS  section  (Figure  3).  The 
DORADO  section  also  taken  on  1 3  August  indicates  that  the 
BL  maximum  is  located  to  the  south  of  the  mooring  location 
(Figure  4).  This  deep  offshore  BL  maximum  disappeared  on 
the  next  day  in  the  DORADO  section  going  through  the  Ml 
mooring  (Figure  4).  The  REMUS  survey  on  15  August 


(Figure  3)  shows  a  shallow  and  weak  BL  signal  coinciding 
with  high  chlorophyll  and  backscatter  signals  around  the  Ml 
mooring,  which  suggests  that  planktonic  dinoflagellates,  not 
zooplankton,  might  be  the  souree  of  this  weak  shallower 
offshore  BL  signal  on  15  August  (Figure  3).  This  might 
indicate  that  with  the  development  of  upwelling,  offshore 
water  masses  with  a  subsurface  layer  of  bioluminescent 
zooplankton  (observed  around  the  Ml  mooring  location) 
were  advected  southward  and  replaced  with  water  masses 
showing  relatively  high  values  of  chlorophyll  fluorescence 
and  backscatter  on  15  August  In  accordance  with  HF  radar 
and  mooring  currents  (Figure  1)  the  presence  of  phytoplank¬ 
ton  in  the  center  of  the  entrance  to  the  bay  might  be  a  result  of 
phytoplankton  advection  from  the  northern  coast  of  the  bay,  due 
to  interaction  of  the  cyclonic  eddy  in  the  bay  and  the  strong 
southward  flow  along  the  entrance  to  the  Bay.  Because  there  is  a 
weak  BL  signal  (Figure  3),  mostly  nonbioluminescent  phyto¬ 
plankton  was  advected  from  the  north  to  the  mooring  Ml 
location.  Surface  predictions  for  10  and  15  August  from  the 
biochemical,  physical  NCOM  ICON  model  indicate  the  deve¬ 
lopment  of  phytoplankton  bloom  in  the  southern  portion  of  the 
bay,  which  extends  from  Moss  Landing  to  the  south  (Figure  5). 
This  is  in  a  good  agreement  with  the  satellite-observed  phyto¬ 
plankton  bloom  on  1 5  August  (Figure  2).  The  model  indicates 
that  one  of  the  reasons  for  this  phytoplankton  bloom  is  an 
increase  of  nutrients  concentrations  (nitrate  and  silicate)  in  the 
southern  portion  of  the  bay  (Figure  6),  where  more  nutrient-rich 
water  is  brought  by  the  southward  flow  and  cyclonic  circulation 
inside  the  bay.  This  increase  in  nutrients  together  with  the 
sluggish  circulation  in  the  southern  portion  of  the  bay  (Figure  2) 
provides  favorable  conditions  for  the  growth  of  diatoms  and 
small  phytoplankton  during  the  development  of  upwelling. 

[20]  Between  10  and  15  August,  the  NCOM  ICON  model 
results  show  the  development  of  a  wide,  well-defined  frontal 
structure  associated  with  the  southward  flow  along  the 
entrance  to  the  bay  (Figure  5).  This  frontal  structure  has  rel¬ 
atively  low  surface  concentrations  of  diatoms  and  zooplank¬ 
ton  (in  comparison  to  the  bay  waters)  and  a  high  concentration 
of  small  phytoplankton  cells  in  comparison  to  adjacent  bay’s 
water  masses  to  the  east  and  offshore  waters  masses  to  the 
west  (Figure  5).  Subsurface  model  results  plotted  along 
REMUS  section  (Figure  7)  also  indicate  a  reduction  in  sub¬ 
surface  diatoms  and  zooplankton  populations  in  the  center  of 
the  entrance  to  the  Bay  (mooring  Ml  area)  as  the  upwelling 
develops.  These  model  predictions  are  in  good  agreement 
with  the  analysis  described  above  of  REMUS  and  DORADO 
observations  during  this  period.  Observations  also  show  a 
reduction  of  zooplankton  and  the  presence  of  small  phyto¬ 
plankton  cells  around  the  Ml  mooring  location  during  the 
upwelling  development  between  1 0  and  1 5  August.  As  it  was 
discussed  above,  high  values  of  chlorophyll  and  backscatter 
in  AUV  observations  at  the  Ml  mooring  location  are  prob¬ 
ably  a  result  of  advection  of  small  phytoplankton  from  the 
northern  part  of  the  bay.  Model  results  support  this:  between 
10  and  15  August,  there  is  a  reduction  in  the  model  small 
phytoplankton  concentration  in  the  northern  area  of  the  bay, 
and,  at  the  same  time,  there  is  an  increase  in  concentration 
along  the  entrance  to  the  bay  (Figure  5).  To  further  verify 
these  proposed  dynamics  during  the  upwelling,  we  present 
derivations  of  an  adjoint  model  to  the  passive  tracer  model. 

[21]  As  shown  in  Appendix  A,  the  adjoint  to  the  passive 
tracer  model  (using  velocities  and  diffusivities  from  the  NCOM 
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Figure  2.  (a)  HF  radar  surface  currents  averaged  over  3  days  of  upwelling  (15-17  August  2003). 
(b)  ADCP  observed  subsurface  profiles  of  the  velocity  components  at  the  Ml  mooring.  Profiles  arc  aver¬ 
aged  over  3  days  of  upwelling.  U  is  the  eastward  component  of  velocity,  and  V  is  the  northward  compo¬ 
nent  of  velocity,  (c)  MODIS-AQUA  sea  surface  temperature  (SST)  and  (d)  surface  chlorophyll  on  15 
August  2003. 


ICON  model)  shows  where  the  model  water  masses  originate 
before  being  circulated  to  the  area  of  interest  [see  also 
Fukumori  el  al,  2004;  Shulman  el  al,  2010].  In  our  ease,  the 
area  of  interest  is  the  area  around  the  Ml  mooring  located  in 
the  center  of  the  entrance  to  the  bay.  Figure  8  shows  the 
adjoint  tracer  distributions,  which  are  estimated  for  volume  V 
in  (A2),  which  consists  of  3  by  3  horizontal  grids  (approxi¬ 
mately  area  of  4  km  by  4  km)  around  the  mooring  Ml  down 
to  a  depth  of  25  m  (the  depth  to  which  high  chlorophyll 
and  baekseatter  values  were  observed  on  15  August,  see 
Figure  3).  Figure  8  shows  vertically  integrated  adjoint  tracer 
maps  at  time  t  equal  to  0000  UT  1 5  August.  The  adjoint  tracer 
maps  are  shown  for  48,  24,  and  1 1  h  prior  to  the  0000  UT  1 5 
August.  In  this  case,  the  adjoint  passive  tracer  distributions 
show  areas  from  which  the  model  tracer-tagged  water  mas¬ 
ses,  48,  24  and  1 1  h  prior  to  15  August,  were  advccted  and 
mixed  into  the  target  area  (around  mooring  Ml).  They  show 
that  these  model  water  masses  mostly  originated  from  the 
northern  part  of  the  bay  around  48  and  24  h  prior  to  15 


August,  and  then  mixed  with  the  water  masses  of  the 
southward  flow  along  the  entrance  to  the  bay.  This  supports 
the  conclusions  above  regarding  phytoplankton  and  zoo¬ 
plankton  dynamics  deduced  from  observations  and  the  bio¬ 
chemical  model  predictions. 

4.  Bioluminescence  and  the  Nighttime  Water- 
Leaving  Radiance  Modeling 

[22]  Moline  et  al  [2007]  proposed  an  approach  for  esti¬ 
mating  nighttime  water-leaving  radiance  (BLw)  due  to  BL 
stimulation  at  depth  which  used  measured  Inherent  Optical 
Properties  (lOPs)  and  BL  to  propagate  bioluminescence 
light  from  depth  to  surface.  However,  the  Moline  et  al. 
[2007]  methodology  does  not  provide  the  capability  to 
forecast  the  nighttime  water-leaving  radiance  on  time  scales 
1-5  days  (only  under  assumption  of  persistence  of  observed 
conditions).  Forecasting  models  of  bioluminescence  inten¬ 
sity  (as  for  example  proposed  by  Shulman  et  al.  [2003, 
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Figure  3.  AUV  REMUS  observed  chlorophyll,  backscattering,  and  bioluminescence  during  11-15 
August.  Solid  vertical  lines  indicate  location  of  the  Ml  mooring. 


2005]),  as  well  as  forecasting  models  of  lOPs  properties  of 
the  water  column  are  needed. 

[23]  In  the  present  study,  the  BL  model  described  in 
section  2.2  is  initialized  on  14  August  by  using  BL  data 


from  four  AUVs  sections:  DORADO  sections  taken  13  and 
14  August,  and  AUV  REMUS  sections  taken  on  14  August 
(Figure  1).  The  initialization  procedure  is  described  in  detail 
by  Shulman  et  al  [2003,  2005].  After  initialization,  the  BL 


Figure  4.  AUV  DORADO  observed  chlorophyll,  backscattering,  and  bioluminescence  on  13  and 
14  August. 
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Figure  5.  Surface  predictions  for  10  and  15  August  from  the  coupled  biochemical,  physical  NCOM 
ICON  model. 


dynamics  are  predicted  forward  in  time  by  using  the 
advection-diffijsion-reaction  model.  Figure  9  shows  surface 
and  subsurface  predictions  of  BL  intensity  on  15  August 
(I  day  later  after  the  model  initialization,  values  are  nor¬ 
malized  by  10g  photons/s).  The  predicted  BL  distributions  in 
the  Bay  are  in  good  agreement  with  the  above  analysis  of 
BL  observations  and  the  biochemical,  physical  NCOM 
ICON  model  predictions.  There  are  high  levels  of  inshore 
modeled  BL  (Figure  9),  which  correlate  with  the  observed 
BL  maxima  and  high  concentrations  of  chlorophyll.  At  the 
same  time,  the  biggest  differences  with  observations  are  in 
BL  distributions  along  the  entrance  to  the  bay  and  offshore. 
The  model  BL  shows  high  values  along  the  entrance  to  the 


bay  (along  the  frontal  structure  discussed  in  the  biochemical 
model  predictions).  However,  as  it  was  discussed  above, 
observations  show  a  very  weak  BL  signal  in  the  area  around 
the  Ml  mooring  location  on  15  August.  It  is  clear  that  the 
high  value  of  model  BL  are  a  result  of  the  advection  by  the 
southward  flow  of  BL  intensity  from  the  northern  part  of 
the  bay.  However,  this  advection  of  bioluminescent  phyto¬ 
plankton  did  not  happen  in  reality  as  illustrated  in  section  3. 
Why  were  bioluminescent  dinoflagellates  not  advected  by 
the  southward  flow?  In  accordance  with  the  biochemical 
model  (Figure  6),  nutrients  were  in  abundance  along  the 
entrance  to  the  bay  (carried  by  the  southward  flow).  At  the 
same  time,  the  entrance  to  the  bay  has  very  strong  currents, 
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which  create  unfavorable  conditions  for  growth  and  survival 
of  dino flagellates  [Jones  and  Gowen,  1990].  This  suggests 
that  advective  and  diffusive  processes  alone  cannot  explain 
many  observed  features  of  spatial  and  temporal  variability  of 
the  BL  intensity.  Modeling  of  behavioral  dynamics  of  bio- 
luminescent  organisms,  as  well  as  modeling  of  sources  and 
sink  terms  representing  ecological  interactions  controlling 
the  bioluminescence  [Benoit-Bird  et  al .,  2010;  Jones  and 
Go  wen,  1990],  should  be  added  to  BL  modeling  method¬ 
ology  in  the  future. 

[24]  Dynamical,  predictive  biochemical  and  biolumines¬ 
cence  intensity  models  provide  the  possibility  to  model  and 
forecast  the  nighttime  water-leaving  radiances  due  to  stim¬ 
ulation  of  BL  at  depth.  We  used  the  constituents  from 
the  NCOM  ICON  biochemical  model  to  estimate  lOPs 
(absorption  and  backscattering)  based  on  the  methodology 
outlined  by  Fujii  et  al.  [2007]. 

[25]  If  the  intensity  of  the  light  from  the  stimulated  BL 
and  lOPs  are  known,  the  propagation  of  the  light  to  the 
surface  can  be  estimated  with  the  radiative  transfer  models 
like  Hydrolight  or  a  reduced  version  of  Hydrolight-Ecolight 
model  [Mobley  and  Sundman ,  2001a,  2001b].  However, 
in  both  cases,  the  use  of  these  models  with  the  coupled, 
biochemical,  physical,  nested,  data  assimilative  models  are 
computationally  expensive.  In  the  present  study  we  esti¬ 
mated  the  propagation  of  light  from  the  BL  source  to  the 
surface  by  inverting  the  Penta  et  al.  [2008]  scheme  which  is 
used  in  the  biochemical  model  for  attenuating  the  PAR  with 
the  depth  (see  section  2.2). 

[26]  Figure  10  (right)  shows  water-leaving  radiance 
(BLw)  at  the  surface  due  to  stimulation  of  the  modeled  BL 
intensity  over  the  entire  model  domain  at  different  depths 
(5,  15,  and  25  m)  on  15  August,  The  modeled  BL  3-D  dis¬ 
tributions  (normalized  by  1CT)  are  shown  in  Figure  10  (mid¬ 
dle)  for  different  depths  of  stimulations  (see  also  Figure  9  and 
above  discussions  about  the  BL  model  results  on  15  August). 
Figure  10  (left)  shows  a  sum  of  a  (absorption)  and  bf,  (back- 


scattering)  averaged  from  the  depth  of  BL  stimulation  to  the 
surface.  Values  of  a  and  bf,  are  estimated  from  the  biochemical 
model  in  accordance  with  Fujii  et  al.  [2007]. 

[27]  There  are  high  values  of  estimated  water-leaving 
radiance  in  the  areas  close  to  the  coast  (Figure  10),  which 
correlates  with  observed  and  model  predicted  BL  maxima 
along  the  coastline.  At  the  same  time,  Figure  10  also  shows 
high  values  of  water-leaving  radiances  along  the  entrance  to 
the  bay  (along  the  frontal  structure  discussed  above),  which 
is  probably  an  artificial  feature  due  to  above  discussed  dif¬ 
ferences  between  modeled  and  observed  BL  distributions 
along  the  entrance  to  the  bay.  The  model-predicted  artifi¬ 
cially  strong  BL  signal  is  illuminated  in  the  surface  and 
subsurface  by  the  relatively  lower  values  of  IOPs  along  the 
entrance  to  the  bay  (Figure  10).  These  low  IOP  values  are 
associated  with  the  development  of  the  frontal  structure  with 
relatively  (to  the  bay  and  offshore)  clear  water  masses  of 
the  southward  flowing  jet  (Figures  2  and  5). 

5.  Conclusions 

[28]  Observations  of  physical,  bio-optical  properties  (includ¬ 
ing  bioluminescence)  together  with  results  from  dynamical 
biochemical,  physical  (NCOM  ICON)  and  bioluminescence 
models  are  used  to  interpret  the  development  of  the  upwelling 
event  during  August  2003  AOSN  11  field  experiment. 

[29]  Satellite  observations  and  the  NCOM  ICON  model 
show  the  development  of  the  phytoplankton  bloom  in  the 
southern  portion  of  Monterey  Bay.  The  model  results  sug¬ 
gest  that  one  of  the  reasons  for  development  of  this  phyto¬ 
plankton  bloom  is  the  increase  of  nutrients  concentrations 
(nitrate  and  silicate)  in  the  southern  portion  of  the  Bay 
where  nutrient-rich  water  masses  are  brought  in  by  the 
southward  flow  and  cyclonic  circulation  inside  the  Bay. 
This  increase  in  nutrients  together  with  the  sluggish  circu¬ 
lation  in  the  southern  portion  of  the  Bay  provides  favorable 
conditions  for  the  phytoplankton  growth. 
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Figure  7.  Subsurface  model  predictions  along  the  REMUS  transect  on  10  and  15  August.  Solid  vertical 
lines  indicate  location  of  the  MI  mooring. 


[30]  Inshore  AUV  observations  show  consistent  coinci¬ 
dence  of  chlorophyll,  backscatter,  and  bioluminescence 
maxima  during  upwelling  development.  Offshore  AUV 
observations  (taken  at  the  entrance  to  the  Bay)  show  deeper 
bioluminescence  maxima  below  the  surface  layers  of  high 
chlorophyll  and  backscatter  values  during  the  earlier  stage 
of  the  upwelling  development.  These  observations  lead  to 
the  conclusion  that  inshore  bioluminescence  maxima  are 
associated  with  the  phytoplankton  (dinoflagellates),  while 
offshore  BL  maxima  are  due  to  larger  zooplankton  which  is 
in  agreement  with  general  differences  in  flash  kinetics 
between  planktonic  dinoflagellates  and  zooplankton  pre¬ 
sented  by  Moline  et  al  [2009].  The  observed  deep  offshore 
BL  maximum  disappeared  during  the  upwelling  develop¬ 
ment  and  became  a  shallower  and  much  weaker  signal 
coinciding  with  high  chlorophyll  and  backscatter  values 


offshore.  Observations  together  with  results  from  the 
NCOM  ICON  model  and  distributions  of  the  passive  tracer 
adjoint  model  suggest  that  with  the  development  of 
upwelling,  the  offshore  water  masses  with  the  subsurface 
layer  of  bioluminescent  zooplankton  were  advected  south¬ 
ward  and  replaced  with  water  masses  showing  relatively 
high  values  of  chlorophyll  fluorescence  and  backscatter. 
This  high  presence  of  the  phytoplankton  at  the  entrance  to 
the  bay  is  a  result  of  its  advection  from  the  northern  coast  of 
the  bay  by  the  strong  southward  flow.  Because  there  is  a 
weak  observed  BL  signal,  mostly  nonbioluminescent  phy¬ 
toplankton  was  advected  from  the  north. 

[31]  Results  from  the  dynamical  bioluminescence  model 
show  agreement  with  the  observed  inshore  BL  observations. 
At  the  same  time,  the  model  shows  high  values  of  BL 
intensity  along  the  entrance  to  the  bay  (along  the  path  of  the 
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Figure  8.  Adjoint  passive  tracer  distributions  on  15  August. 


southward  flow);  while  observations  show  that  biolumi- 
nescent  dinoflagellates  mostly  escaped  the  southward 
advection  from  the  northern  part  of  the  bay.  This  indicates 
that  the  BL  modeling  approach  (which  is  based  on  advcctive 
and  diffusive  processes)  should  include  the  modeling  of 
behavioral  dynamics  of  bioluminescent  organisms,  as  well 
as  the  modeling  of  sources  and  sink  terms,  representing 
ecological  interactions  controlling  the  biolumincscencc. 


[32]  The  biochemical  and  bio  luminescence  models  pre¬ 
sented  in  this  study  are  combined  into  a  methodology  for 
estimating  the  nighttime  water-leaving  radiance  due  to  stim¬ 
ulated  bioluminescence.  The  results  show  high  values  of 
estimated  water-leaving  radiance  in  areas  along  the  coastline, 
where  high  values  of  BL  were  observed  and  model  predicted. 
At  the  same  time,  estimated  water-leaving  radiances  have 
high  offshore  values  (along  the  entrance  to  the  Bay),  which 


A  MODEL  BIOLUMINESCENCE 


TIME  (HOUR) 

Figure  9.  (a)  Surface  and  (b)  subsurface  (along  the  REMUS  transect)  model -predicted  BL  distributions. 
Solid  vertical  line  indicates  location  of  the  M 1  mooring. 
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Figure  10.  (right)  Water-leaving  radiance  at  the  surface  due  to  stimulation  of  the  modeled  BL  intensity 
at  different  depths  (5,  15,  and  2  5m)  on  15  August,  (middle)  The  modeled  BL  intensity  for  different 
depths  of  stimulations,  and  (left)  a  sum  of  a  (absorption)  and  bh  (backscattering)  averaged  from  the  depth 
of  BL  stimulation  to  the  surface. 


are  comparable  to  the  coastal  distributions.  These  off¬ 
shore  high  values  of  water-leaving  radiance  are  a  result  of 
artificial,  model-predicted  BL  maximum  along  the  entrance 
to  the  Bay. 

Appendix  A:  Passive  Tracer  and  Its  Adjoint 

[33]  Consider  the  passive  tracer  equation  for  concentration 
C(x,  y,  z,  t) 


With  initial  conditions  at  time  t  =  to 

C  =  C0, 

where  diffusivities  (A  and  K.)  and  velocities  (u,v,w)  are  from 
the  NCOM  ICON  model  described  in  section  2.2. 

[34]  Let  us  consider  the  following  objective  function  J  at 
time  t  >  to 


(A2) 
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where  V  is  a  particular  subdomain  (target  area)  of  the 
modeling  domain,  r  is  the  location  in  the  model  domain 
with  coordinates  (x,y,z),  and  dr  is  a  volume  element. 
Therefore,  function  J  is  the  normalized  content  of  tracer  C  in 
the  domain  V  at  time  t. 

[35]  By  using  the  adjoint  for  the  tracer  equation  (Al),  the 
gradient  of  the  function  J  (equation  (A2))  at  time  t  with  respect 
to  the  initial  concentration  Q>  at  time  to,  can  be  estimated 


where  s  is  the  sensitivity,  is  the  gradient  of  J  (at  time  t)  with 
respect  to  initial  conditions  ^0.  Sensitivity,  s,  is  a  function  of  x, 
y,  z  and  times,  to  and  t,  and  can  be  estimated  by  seeding  the 
adjoint  variable  with  a  unit  value  at  each  grid  point  in  the 
volume  V  at  time  t,  and  integrating  the  adjoint  of  the  tracer 
model  backward  in  time  to  time  to  [Fukumori  et  al. ,  2004; 
Shulman  et  al l,  2010].  The  function  s(x,  y,  z,  t)  is  called  the 
adjoint  tracer  sensitivity,  as  well  as  the  adjoint  tracer  distri¬ 
bution  (because  the  function  s,  is  the  result  of  the  adjoint  tracer 
model  integration). 

[36]  Lets  introduce  some  finite  perturbation  AC0  at  loca¬ 
tion  Xo  =  {xo,  yo,  Zq}  to  the  initial  concentration  C0  at  time 
to,  according  to  (A3)  we  would  have 

AJ  =  s(xo,  y0, zo» to,  t)  •  ACft.  (A4) 

According  to  (A2)  and  (A4),  the  adjoint  tracer  distribution 
s(xo,  yo,  Zo,  to,  t)  represents  a  fraction  of  tracer  AC0,  which 
makes  its  way  to  the  volume  V  from  time  to  to  time  t.  Due  to 
the  linearity  of  the  passive  tracer  and  its  adjoint  problems, 
the  adjoint  tracer  distribution  s(x,  y,  z,  to,  t)  will  represent 
the  fraction  of  the  tracer-tagged  water  that  makes  its  way 
from  location  (x,  y,  z)  at  time  to  to  the  volume  V  at  time  t 
[see  also  Fukumori  et  al .,  2004].  Therefore,  the  adjoint 
tracer  distribution  provides  information  on  the  model  tracer 
history  and  identifies  origin  and  pathways  of  the  model 
tracer-tagged  water  masses  in  the  past,  which  circulated  into 
the  target  area,  and  therefore  contain  useful  information 
about  model  circulation  patterns  and  the  propagation  of 
information  within  the  system. 
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